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A  TWO-DIMENSIONAL  COMPRESSIBLE,  NONLINEAR 
STABILITY  ANALYSIS  WITH  APPLICATIONS  TO 
THE  STUDY  OF  BOUNDARY  LAYER  TRANSITION 

I.  INTRODUCTION 

Recent  analytical  efforts  to  gain  a  more  comprehensive  understanding 
of  boundary  layer  stability  have  centered  on  (1)  the  coupling  between 
the  principle  linear  wave  (the  Tollmien-Schlichting  wave)  and  its  principle 
harmonics  (for  example,  [1])  and  (2)  on  the  development  of  computer  codes 
to  study  the  nonlinear  effects  of  disturbances  propagating  through  an 
incompressible  boundary  layer  field  (for  example.  Refs. 2,  3,  4).  Such 
investigations  extended  the  earlier  incompressible  and  compressible  linear 
stability  analyses  (for  example.  Refs.  5-9)  which  analyzed  the  linear 
stability  problem  both  from  spacial  and  temporal  growth  processes. 

The  recent  efforts  have  been  motivated  by  analysis10  and  experiments11 »^»13 
which  indicated  that  the  process  of  disturbance  growth  and  instability  leading 
to  transition  could  be  subdivided  into  distinct  regions, where  first  linear, 
and  then  nonlinear  and  three  dimensional  effects  predominate. 

Two  dimensional  nonlinear  codes  can  thus  be  extremely  important  tools 
If  one  is  to  properly  stimulate  the  three  dimensional  bursting  phenomenon 
of  transition  to  turbulence,  since  such  codes  can  provide  the  initial  condi¬ 
tions  for  such  analysis  (simitar  to  the  effect  of  the  linear  codes  which 
provide  the  initial  conditions  for  the  nonlinear  analyses). 

While  the  nonlinear  codes  available  represent  a  significant  advance 
in  computational  capability,  they  have  thus  far  been  restricted  to  incompres¬ 
sible  flow  calculations.  Such  programs  are  unable  to  analyze  the  effects  of 


Mach  number,  wall  temperature,  imposed  disturbance  propagation  velocity,  etc. 

Experiments14,  ^  have  indicated  that  imposed  acoustical  disturbances 
effect  transition  in  a  low  speed  flow  in  a  manner  similar  to  free  stream 
turbulence,  even  though  the  propagation  velocity  of  the  disturbances  are 
as  much  as  two  orders  of  magnitude  greater.  The  frequency  range  which 
effects  transition  at  a  given  Reynolds  number  is  similar  to  the  range 
associated  with  freestream  turbulence,  even  though  there  is  a  mismatch 
in  freestream  propagation  velocity  (and  thus  in  freestream  wave  number). 

The  purpose  of  the  study  at  New  York  University  for  the  past  three  years 
was  to  ascertain  the  mechanisms  which  cause  acoustical  disturbances  to  effect 
transition  in  this  manner  and  to  develop  a  computer  code  which  includes  both 
nonlinear  effects  and  the  effects  of  compressibility.  The  unsteady,  two- 
dimensional,  compressible  second  order  boundary  layer  equations  (retaining 
terms  of  the  order  of  the  reciprocal  of  Reynolds  number  squared)  have 
been  computed  numerically1®.  The  system,  consistent  with  the  Navier  Stokes 
equations,  have  been  solved  by  MacCormack  explicit  scheme  with  time 
splitting17.  In  Section  II  the  equations  utilized  and  the  method  of  analy¬ 
sis  is  developed.  Results  utilizing  the  program  are  presented  in  Section  III. 
A  description  of  the  program's  capabilities  can  be  found  in  Section  IV. 


II.  METHOD  OF  ANALYSIS 

A.  Equations  of  Motion 

The  equations  of  motion,  the  unsteady,  compressible  Navier  Stokes 
equations,  were  written  in  terms  of  the  nondimensional  parameters  consistent 
with  the  boundary  layer  equations  (but  terms  of  order  of  the  reciprocal 
of  the  Reynolds  number  squared  were  retained).  The  equations  were  trans¬ 
formed  by  a  stretching  parameter  n  which  packs  most  points  in  the  boundary 
layer  but  allows  for  a  significant  distribution  of  data  in  the  freestream. 
The  stretching  utilized  was; 

n  *  1  -  exp  (  -ay/6*) 


where  3*  =  ^  ,  y  *  ^-jRe^,  and  6o*(x)  is  the  initial  boundary 

layer  thickness  distribution  (the  distribution  prior  to  the  initiation 
of  perturbed  flow).  ReL  is  the  freestream  Reynolds  number  based  on 
characteristic  length  L  (for  example,  plate  length). 

The  resulting  equation  systeu^,  are  rewritten  here: 

■^  +  -|x(pi3)+ln  (pG)  lr+(1"n)  fj  !?r“  "  0  (1) 


3_ 

at 


o 


3 

3n 


(pu)  +  -|y  (55“  +  p)+fj-  (pu®  +  pj  (1  -n) 

[5  0  -«)  |jj-l  ♦  ♦. 

fe-  <5?>  A  IT  ♦  in  (5=5)  ft-  *  ('  -nlf.'ln  (W*  +  *  ■*> 
—  pzd  -  n)  5^  [5(1-  n)  ji]  -  |-(i  -  n)  |f  If  [5  |jj-  +  5  §  |f  ] 

*  fe®.  0 -ft  tflBO-nJfr 

ft  (*)’  ♦  L  <5=">  ♦  f-  (55">  |_(55")  • 


-  tt2  (1  -  n)  3_£y(l  -  n)  3H  1 

£*»  3n  3n 


(2) 

(3) 


(4) 


+ 


where 


x  =  x.  y  *  y  •  v'Rij ;  t  =  t  u 
L’  '  L  L  t“ 

5*ib  5  *i:  *  p 


P _ .  -  .  p 

P  iF  *  p  P 
r®  •  r® 


K 

12 


¥  ♦ 

o» 


E  _  CvT  .  u2  +  v2 

iF  7r'+  ~~2u2 
®  ®  0» 


V. 


T«  +  198.6  j  .  RT  Re,  ■  pmiy. 

T  +  198.6  *  T  uj  L 


_  1 _  ,  4  3  3u  £  a  3 

■Re[”  1  3  3x  3x  "  3  j*  3x 


ffr  IsrCu  (l-n)-g  ] 

T"  lx  In  (p  f£}  ’  I  ^  "Uln"^  0  *  n)  -"In1 

0  "«>  fr  I) 


4. 

3 


P 


'  n'  y  “5n  a?  J 

3^  (-  _5n  _iy.\  +  i  3n  i_  /:  Ja 
3x  'p  3X  3n;  3  3x  3n  3x  an' 

I1  -  -0  It  g  ft  i» 


§*  *  ”  3n  ax  an' 

1  r  3_  t~  3v  \  ,  3 

*S[C3X(l1  aT}  +  ax  ai 

♦  fc<& 

< cfr  ♦  -ftlfi  csM 


ft  ft  (5f» 


'A  wn  VA  01 1  ^ 

i,  (  iS  i£)  +  (ill)  i_  (inii  )] 
ax  '  ax  3n;  '3X'  3n  'ax  3n  ;J 

>  an  a _ r-  3H  ;  an  3Hi 

ax  an  J  Ly  3X  a3T  anJ 

ia  _au  2  i  3v  +  in iij? 

ax  3n  *  3  M\  ax  ax  aV 


/3u 

'ax 


2  a(l  -  n)  -  /  3u  _3a  i“  ^  ii 

3  j  *  u  '  ax  +  ax  3n  ;  an 

+  2  alLz.Il)  -  r  ( |v  3j  3v J  afl 

;  *  w  '  5T  ax  3n;  an 


.  l 


i<ft  ♦  ft  ft>  <3F-*-ft-R) 


-4- 


* 


Equations  (1)  -  (4)  have  been  solved  by  a  MacCormack  predictor-corrector 
explicit  scheme  with  time  splitting^).  At  t  *  o,  a  disturbance  field, 
calculated  utilizing  linear  stability  theory  (see  below),  is  initiated  near 
the  leading  edge  of  a  plate.  The  time  dependent  flowfield  is  then  calculated 
numerically  from  the  unsteady,  compressible  equations,  consistent  with  the 
imposed  initial  and  boundary  conditions. 

B.  Initial  and  Boundary  Conditions 

The  investigation  has  centered  on  a  study  of  the  amplification  or 

dampling  of  acoustical  disturbances  propagated  into  a  boundary  layer. 

14 

The  experiments  of  Spangler  and  Wells  who  utilized  an  air-driven, 

rotating  vane  sound  generator  to  create  the  disturbance  without  pro¬ 
ducing  any  appreciable  turbulence,  were  simulated.  Both  the  frequency 

and  the  intensity  of  the  sound  source  had  been  varied  experimentally. 

In  the  experiments,  a  low  velocity  boundary  layer  channel  was  run  at  a 
unit  Reynolds  number  of  2.4  x  1 0s/f t. ,  the  channel  wall  representing 
the  flat  plate.  Measurements  of  transition  occurred  at  distances  on 
the  order  of  10  -  20  feet,  and  thus  the  length  Reynolds  numbers  of 
interest  were  in  the  106  to  107  range. 

In  order  to  model  this  problems,  a  set  of  initial  and  boundary 
conditions  must  be  established,  consistent  with  the  propagation  of 
acoustical  disturbances  and  at  the  same  time  consistent  with  the 
set  of  differential  equations  utilized.  The  initial  data,  consistent 
with  equations  (1)  through  (5)  under  steady  state  conditions,  was  origi¬ 
nally  established  utilizing  numerical  techniques  for  subsonic  boundary 
layer  analysis  with  normal  pressure  gradients  established  by  the 
principal  investigator  and  reported  previously. 16  Thus  the  solution 
t  *  o  was  the  solution  of  the  two  dimensional  steady  boundary  layer 
equations  with  normal  momentum  equation  in^uded. 
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Consistent  with  Refs.  2  and  3  similar  solutions  for  the  initial 
profiles  were  studied.  Since  no  problems  were  encountered,  the 
final  program  now  utilizes  such  solutions  for  the  initial  field. 

At  the  wall  (n  =  o),  for  a  rigid  body,  one  can  neglect  the  effect 
of  the  wave  on  temperature,  and  for  an  adiabatic  wall  one  can  establish 
the  following  relations; 
at  n  =  0  (x  >  x*) 
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The  outer  boundary  condition  (n  -*•  1)  is  established  by  allowing  the 
wave  to  travel  as  a  plane  wave  of  arbitrary  speed  (for  the  acoustic 
wave,  the  speed  of  sound  C^). 

Much  effort  was  expended  on  the  determination  of  a  proper  downstream 
boundary  condition.  It  was  found,  consistent  with  the  results  of 
Fasel2,  that  the  boundary  condition  that  yields  the  least  upstream 
influence,  and  is  thus  superior  to  other  possible  ones  (including  a 
non- reflective  condition)  is: 


32 


u* (xf,  y,  t) 
3x2 


a2u' 


where  a  *  fL/Vph  (f  is  the  frequency  and  Vph  the  phase  velocity),  u' 
is  the  perturbation  quantity  (u  -  uQ).  and  x^  is  the  downstream 
boundary.  This  condition  says  that  at  the  downstream  boundary,  the 
disturbance  has  a  periodic  form.  At  t  a  o,  a  disturbance  is  initiated 
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at  x  ^  o  so  that  the  velocity  field  at  x  *  x*  =*  o  and  t  ■>  o  can 
be  written  as 

u(x*.  y,  t)  =  u0(x*.  y)  +  fi(y)  sin  (ut) 
where  uQ(x*,  y)  is  the  velocity  profile  at  t  *  o,  u  *  (f  being 

the  frequency  of  the  imposed  disturbance,  in  cycles  per  second)  and 
f i (y)  is  the  disturbance  profile  near  the  leading  edge. 

The  form  of  fx  (y)  can  in  general  be  prescribed 
and  then  the  eigenfunction  associated  with  the  eigenvalue  problem 
detailed  below  can  be  superposed  to  reconstruct  fi(y).  The  complete 
determination  of  these  eigenfunctions  is  quite  complicated  in  general 
depending  on  the  initial  boundary  layer  flow,  and  leads  to  a  signifi¬ 
cant  numerical  problem  in  linear  stability  theory.  Such  complications 
as  the  existence  of  a  continuous  spectrum  of  eigenvalues,  the  presence  of 
a  singularity  in  the  linearized  eigenvalue  problem,  and  the  existence  of 
viscous-type  eigenstates  which  behave  singularly  in  the  inviscid  limit. 
[Refs.  18-20]  are  indicative  of  studies  of  such  problems. 

The  other  boundary  conditions  at  x  *  x*  are  consistent  with  linear 
theory  and  are  to  be  determined  once  the  form  of  fx(y)  is  known. 

In  order  to  derive  appropriate  perturbation  profiles  for  use  as 
the  initial  boundary  condition  in  our  program  (at  x  *  x*)  the  compressible 
analogue  of  the  Orr-Sommerfeld  equations  must  be  solved.  For  purposes 
of  the  present  study,  we  have  only  considered  these  equations  in  the 
inviscid  limit  and  under  the  additional  assumption  of  no  temperature 
gradient  in  the  mean  flow. 

Thus, nondimensionali zing  all  lengths  by  the  boundary  layer  thick¬ 
ness  6  at  a  suitable  station  along  the  plate,  all  velocities  by  the 
free  stream  speed  of  sound,  and  the  pressure  by  ,  the  linearized 
perturbation  system  considered  is; 
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where  bars  have  been  dropped  and  subscript  zero  refers  to  the  profile 

at  t  =  o.  The  appropriate  boundary  conditions  are  discussed  below. 

As  is  customary  the  normal  mode  type  decomposition  has  been  used, 

•  •  •  -  •  *  ■  • 

and  hence  we  note  that  p,  u,  V,  p,  T  are  functions  of  y,  the 
distance  normal  to  the  plate.  The  parameters  u  and  k,  the  dimension¬ 
less  frequency  and  wave  number  respectively  are  given  by 
a)  =  £25/  ,  k  *  1(6 

Coo 

where  £2  and  fc  are  dimensional  quantities.  The  mean  flow  quantities 
“o  =  M-uB  and  “i  ■  5.6  (M^S) 

are  the  dimensionless  velocity  and  velocity  gradient  respectively  where 
MB  is  the  freestream  Mach  number,  ug  is  the  dimensionless  (with  respect  to 
U,,)  Blasius  profile  and  the  factor  5.6  results  from  changing  the  normal 
coordinate  from  n  to  y. 

Several  simplifications  can  be  made  in  equations  (6  -TO)*  We  have 
solved  for  p,  u,  v  and  T  in  terms  of  the  pressure  p  and  find  that  p  must 
satisfy  the  second  order  equation 

P"  -  P'  +  [  (“9.  *  -  h>)2  -  k2]  p  «  o 

°  “  T0 
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we  note  that  the  second  condition  is  a  normalization  condition. 

Equations  (11)  -  (12)  now  provide  an  eigenvalue  problem  for 
the  eigenvalue  k.  Once  the  possible  values  of  k  are  determined, 
the  profiles  x  *  x*  (i.e.  the  fi(y))  can  be  determined  and  thus  the 
initial  and  boundary  conditions  for  the  calculation  are  entirely 
prescribed.  To  treat  this  eigenvalue  problem  a  program  developed  by 
Mack  Ref.  (23)  has  been  used.  The  program  has  the  capability  of 
treating  cases  where  the  phase  velocity  of  the  disturbance  equals 
the  base  flow  velocity  (critical  layer)  by  deformation  of  the 
integration  contour  into  the  complex  plane. 
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III.  RESULTS  UTILIZING  PROGRAM 


The  nonlinear  computer  program  for  the  solution  of  equations  (1)  -  (4) 
was  first  run  for  cases  corresponding  to  the  Reynolds  number  and  frequency 
range  of  the  experiments  of  Reference  14.  Figure  (1)  presents  the  results 
of  three  such  computations  for  the  forced  response  of  the  boundary  layer 
to  acoustic  waves  of  different  frequencies  at  a  freestream  R.M.S.  level  of 
.3#.  The  results  are  in  line  with  Tollmien  Schlicting  amplication  rates. 

An  investigation  was  initiated  to  determine  what  the  propagation 
speed  inside  the  boundary  layer  was.  In  all  cases,  while  the  disturbances 
were  propagating  with  the  speed  of  sound  along  the  freestream,  the  waves 
within  the  boundary  layer  were  propagating  at  a  speed  of  the  order  of 
the  freestream  speed,  and  thus  indeed,  the  waves  are  the  classical  Tollmien 
Schlichting  waves.  Such  a  result  is  presented  in  Figure  2. 

What  is  essentially  occurring  is  that  the  wave  propagating  along  the 
outer  edge  has  little  effect  on  the  boundary  layer  development.  Instead, 
the  major  effect  is  the  profile  at  the  leading  edge,  and  thus,  the  effect 
of  the  acoustic  wave  is  only  to  set  up  the  initial  disturbance  field  (i.e. , 
near  x  *  o),  as  described  previously.  The  effect  of  acoustic  waves  on 
transition  is,  thus  similar  to  the  effect  of  other  imposed  perturbations 
such  as  freestream  turbulence. 

Another  numerical  experiment  was  initiated  to  test  this  result.  The 
program  has  the  capability  to  adnlt  freestream  disturbances  of  any  speed. 
Waves  were  propagated  at  the  freestream  velocity,  instead  of  at  the 
acoustic  speed,  and  all  other  parameters  were  held  constant.  The  dashed 
line  in  Figure  2  indicates  that  the  effect  on  the  boundary  layer  (y/6  <  1) 
is  insignificant,  and  the  outer  boundary  condition  only  effects  the  flow- 
field  through  its  amplitude. 
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The  boundary  layer  response  to  an  imposed  disturbance  is  thus  deter¬ 
mined  completely  once  the  initial  profile  (from  linear  stability  theory), 
is  found  through  a  complete  eigen-value  search.  Nonlinear  programs  should 
therefore  be  coupled  to  linear  stability  programs  which  determine  the 
linear  (upstream  boundary  condition)  profiles  deduced  from  the  multitude 
of  eigen  values  which  have  been  found  to  exist. 
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IV.  PROGRAM  CAPABILITIES 


A  computer  program  which  can  investigate  the  nonlinear  response  of  a 
compressible  boundary  layer  on  a  flat  plate  to  any  imposed  external  distur¬ 
bance  has  been  produced.  The  program  first  computes  the  boundary  condition 
near  the  leading  edge  utilizing  linear  stability  theory  and  then  determines 
the  disturbance  field  downstream  as  a  function  of  time,  utilizing  the 
MacCormack  time  splitting  explicit  technique. 

The  program  originally  created  to  provide  a  fundamental  understanding 
of  boundary  layer  response  to  acoustical  disturbances  has  been  modified  so 
that  a  wide  range  of  parameters  can  be  varied.  The  effects  of  Mach  number 
and  wall  conditions  can  be  computed  utilizing  the  program. 
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FIGURE  2. 


PHASE  SPEED  OF  WAVES  INSIDE  BOUNDARY  LAYER  UNDER  THE  INFLUENCE 
OF  ACOUSTICAL  DISTURBANCE  (DASHED  LINE  INDICATES  PHASE  SPEED 
WHEN  FREESTREAM  WAVE  PROPAGATES  AT  U  ) 
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